Analytical solution of a one-dimensional Ising model with zero temperature dynamics 
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The one-dimensional Ising model with nearest neighbour interactions and the zero-temperature 
\ dynamics recently considered by Lefevre and Dean -J. Phys. A: Math. Gen. 34, L213 (2001)- is 

. investigated. By introducing a particle-hole description, in which the holes are associated to the 

domain walls of the Ising model, an analytical solution is obtained. The result for the asymptotic 
energy agrees with that found in the mean field approximation. 
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I. INTRODUCTION 



o 

^ , Ising systems are often used as simple models in many different fields of statistical physics. In particular, the 
one-dimensional Ising model with nearest neighbour interactions and Glauber dynamics pi is ubiquitous. It shows 
glass-like behaviour, including non-exponential relaxation a laboratory glass transition in cooling processes [^||, 

^ • hysteresis effects in thermal cycling experiments [Q, and aging effects at low temperatures |l^-|l^. When a weak 
c/2 \ oscillating external field is applied, the amplitude of the induced magnetization presents a maximum as a function of 
the temperature, a behaviour resembling stochastic resonance p^ , p^ . 

Other kind of important Ising systems are those with "facilitated" dynamics. They were introduced in a pioneering 



work by Fredrickson and Andersen p6 17|, in the context of the relaxation of structural glasses. In these models, 
a given spin can only flip if its nearest neighbours are in a certain subset of all their possible configurations. The 
^ ] one-dimensional Ising model with facilitated dynamics (ISFM) has been extensively studied p^-p^. Although, in 
Q ■ general, it is not possible to find an analytical solution of these facilitated Ising models, the ISFM has been exactly 
O ' solved at zero temperature [^^. Very recently, the ISFM has been applied to study the compaction of vibrated 
granular systems. A hole (particle) is associated to any spin in its excited (ground) state. Trying to mimic what is 
I done in real experiments |2|,^, the tapping process is simulated in the following way. First, the system freely relaxes 
^ : to a metastable state (n = 0) following the dynamics at zero temperature, in which the processes decreasing the 
^sD ' density of particles are forbidden and any configuration with all the holes being isolated is an absorbent state of the 
\ dynamics ||2^]. Then, the system is tapped, allowing the density to decrease. Afterwards, a new free relaxation at zero 
temperature is done, reaching another metastable state [n = 1). By repeating this process, a chain [n = 0, 1, 2, . . .) 
of metastable configurations is obtained, with the density of the system increasing as a function of n. 

An analogous approach to the problem of granular compaction has been carried out by Lefevre and Dean |p8| , p9| . 
^ ■ Their model is again the one-dimensional Ising system with nearest neighbour interactions. Nevertheless, the usual 
^ Glauber dynamics [Q is not a good choice for simulating tapping processes. This is because there are no metastable 
states, and at zero temperature the system always reaches the perfectly ordered ferromagnetic phase. Then, other 
^ , zero temperature single spin flip dynamics was considered, in which only the elementary events lowering the energy 
I ' are possible. In this way, all the states composed of domains of / parallel spins with length I > 2 are metastable, 
' O i. e., they are absorbent states for this dynamics. In such states, there is no spin antiparallel to both of its nearest 
neighbours, which are the only spins being able to flip. With this falling dynamics, the tapping process is simulated 
in the same way as described above. 

In ref. |p8|| , the authors claimed that the new zero temperature dynamics does not seem amenable to analytic 
solution. In fact, the usual way of solving the Glauber model, by constructing the hierarchy of equations for the 
two-spin moments, does not work. The aim of this paper is to show that it is possible to solve analytically the time 
;J] ' evolution of the model by going to an equivalent particle- hole description of the Ising system. The holes are associated 
I to the domain walls separating arrays of parallel spins (particles). Within this picture, the metastable states are those 
with all the holes being isolated, as in the one-dimensional facilitated Ising model. 
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The paper is organized as follows. In section || the Ising model is introduced, as well as the equivalent particle- hole 
description. Section [II is devoted to the analytical solution of the dynamics. A closed hierarchy of equations is 
derived for the probabilities of finding r+1 consecutive holes. The solution of this hierarchy is obtained by means 
of a generating function method. In particular, the asymptotic density of holes and energy in the metastable state 
are exactly calculated. The last section contains some final remarks. 



II. THE MODEL 



We consider the one-dimensional Ising model with nearest neighbour interactions and periodic boundary conditions. 
The hamiltonian of the system is 



N 

H = -j'^aia^+i , (1) 



where J > is the coupling constant, having dimensions of energy, N is the number of spins, and (jj = ±1. The time 
evolution of the system is described by a Markov process with single spin flip dynamics. Then, the probability p{(T, t) 
of finding the system in configuration cr = {cTi} at time t obeys the master equation 

^p((T, t)^Yl [W{(T\R,(T)piR,rT, t) - W(RM(^)p{c7, t)] . (2) 

i 

Here Ricr is the configuration obtained from cr by flipping the i-th spin, and W{Ri<j\(7) is the transition rate for 
that process. Following Lefevre and Dean psf , we consider a zero temperature dynamics in which only the spin flips 
lowering the energy are permitted. Namely, we take 

W{R^Ct\(t) - I (1 - <^^-l<T^) (1 - (T^(J^+l) , (3) 

i. e., the transition rate equals a if the i-ih spin is antiparallel to both of its nearest neighbours, and vanishes 
otherwise. While an exact solution of the usual Glauber dynamics can be found for all temperatures from the 
hierarchy of equations for the two-spin moments the same procedure applied to this "falling" dynamics leads to 
a non-closed set of equations, since more complex moments, involving three and four spins, appear in them. 

For the above reason, it is convenient to go to an equivalent particle-hole description of the Ising model. For each 
site i, we define a new variable 

m.= ^-^. (4) 

When rrii — 1 we will say that site i is occupied by a hole, while if to^ = we will refer to site i as being occupied by 
a particle. Thus, there is a hole at site i if the spins i and i + 1 are antiparallel, while a particle corresponds to spins 
i and i + 1 being parallel. It follows that holes are associated to the domain walls separating arrays of parallel spins 
(particles). It is important to note that the number of holes is even for any configuration with periodic boundary 
conditions. In terms of the new variables, the hamiltonian of the system reads 

N N 

jY^{2m,-l) = -JN + 2J^m,. (5) 

1=1 i=l 

We define the average dimensionless energy per spin e as 

N 

JN ^ ' N 



where the angular brackets (. . .)< denote average with p(cr, t). 

In the particle-hole description, the elementary events involve two adjacent sites. When the i-th spin flips, both 
TOi„i and rrii change their state. From Eq. (^, the transition rate w(i?i_ii?im|m) for this process is 

uj(Ri-iRim\m) — ami-irrii . (7) 
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As a only determines the arbitrary time scale, we will set a = 1 in the following. The master equation for the 
probability p(m, t) of finding the system in configuration m at time t reads 



^p(m, t) = ^ \^{rn\Ri-xRirri)-p{Ri-\Rim, t) - Lo{Ri-iRim\m)p{m, t)] . (8) 

i 

In the dynamics given by the transition rates in Eq. (Q), only two nearest neighbour holes can turn into two particles. 
Therefore, it is clear that, after a long enough time period, the system will become stuck in a "metastable" state with 
all the holes being isolated, i. e., surrounded by two particles. Of course, the reached metastable state will depend 
on the initial configuration. This behaviour is reminiscent of the one showed by the ISFM at T = 0. In the latter, 
the system also reaches a metastable state with all the holes isolated, the specific final state depending on the initial 
condition ||2^,^. Nevertheless, the model considered here and the ISFM at zero temperature are not equivalent. 
The elementary dynamical events occurring in each of them are different. While a particle can be adsorbed on any 
empty site with at least one nearest neighbour hole in the ISFM, in the present model two particles must be adsorbed 
simultaneously on two adjacent empty sites of the one-dimensional lattice. 

III. ANALYTICAL SOLUTION OF THE DYNAMICS 

In order to analyze the dynamics of the model we focus on the set of moments 

Dr{t) = {mkruk+i . ■ .mk+r)t , (9) 

with r > I. The local density of holes is given by the first term of this hierarchy, 

Do{t) ^ {mk)t , (10) 

and it is directly related with the energy of the system. From Eq. (^), it is 

e = -l + 2Do. (II) 

Note that we are restricting ourselves to spatially homogeneous situations, so these moments Dr do not depend on 
their first site k. From its own definition, the moment D^- gives the probability of finding r + I consecutive holes 
starting from a given arbitrary site of the lattice. Using the master equation (^ with the transition rates (0) one gets 

^^Dr{t) = -2Dr+iit) - rDr{t) , (12) 
for all r > 0. This hierarchy can be solved by introducing the generating function 

OO J. 

from which all the moments Dr{t) are easily obtained, 

Dr{t) ^ \im^^G{x,t) . (14) 

The hierarchy of equations (|l2|) is equivalent to the following first order partial differential equation for the generating 
function: 

dtG{x,t) = -{2 + x)d^,G{x,t), (15) 
which has to be solved with the initial condition 



Go(x)=G(x,0)=^^A.(0). (16) 



r=0 

By using standard techniques it is easily obtained 
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Gix,t) = Go[{2 + x)e-'-2]. (17) 
For large times the solution approaches the limit 

G(a;,oo) =Go(-2), (18) 

and taking into account Eq. (|l^), 

lim Do(0 = Go(-2) , lim A-(<) = , r > 1 . (19) 

t — >00 t — >QO 

This shows that all the holes are isolated in the metastable state reached by the system in the long time limit. The 
probability of finding r + 1 consecutive holes, which is equal to Dr, vanishes for r > 1. On the other hand, the 
asymptotic value of the density of holes Do{oo) does depend on the initial state, being always smaller than its initial 
value. In this sense, it is worth noting that the hierarchy of equations (^2|) has as a stationary solution any constant 
value of Do, and Dr = ior r > 1. 

Now, we are going to particularize the above results for a specific initial condition. Let us consider that the one- 
dimensional Ising model is initially at equilibrium at a certain temperature T. Then, the probability that spins i and 
i + 1 are antiparallel is given by a function of the temperature a{T), 

expf — •^'^ ) 

where ks is Boltzmann's constant. The value a = 1/2 corresponds to T — > oo, the system is at i = in a completely 
random configuration. From the form of the hamiltonian in the particle-hole description, Eq. (|^), it follows that at 
equilibrium the nik variables are statistically independent, and 

Dr{0) = a"+' . (21) 

By using Eq. (^6|), one gets the initial condition for the generating fmiction, 

Go{x) = ae""" , (22) 

and Eq. (|l|) yields 

G{x, i) = a exp {a [(2 + x)e"* - 2] } . (23) 

Now the expressions for the moments Dr{t) can be directly calculated from Eq. (p^, 

Drit) = a^+'' exp [-rt + a (2e"' - 2)] . (24) 

As discussed in the previous paragraph, all the moments Dr with ?' > 1 vanish in the long time limit. Moreover, the 
asymptotic density of holes is 

£'0(00) = ae^^'^ , (25) 



which depends on the initial temperature of the system. Equivalently, from Eq. ( |1 1| ) the specific dimensionless energy 
e(cxo) in the metastable state reads 

£(00) = -1 + 2ae"2a _ (26) 

The mean field calculation presented by Lefevre and Dean |^,^ leads to the same expression (see also [Q). The 
asymptotic energy £(c!o) is maximal for a = 1/2, i. e., when the system starts from the completely random configura- 
tion. This is in contrast with the behaviour found for the ISFM. In the latter, if the initial condition is taken as the 
equilibrium state for a given value of the temperature T, the asymptotic density of holes is a monotonic function of 
the initial density of holes |2^ . 
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IV. FINAL REMARKS 



In this work we have studied the time evolution of the one-dimensional Ising model with the zero temperature 
dynamics recently considered by Lefevre and Dean [2^p9[| . In this dynamics only transitions decreasing the energy 
of the system are allowed. The difference with the usual zero temperature Glauber dynamics |^] is that flips of spins 
parallel to one of its neighbours and antiparallel to the other one are forbidden. In other words, the possibility of a 
diffusive motion of the domain walls is eliminated. 

When trying to compute analytically the evolution of the system, the usual procedure of constructing the hierarchy 
of differential equations for the two spin moments is not useful. We have introduced a particle-hole description of the 
Ising model, in which the holes correspond to the domain walls (or "defects") between arrays of parallel spins. The 
energy of the system is directly related to the density of holes. Within this description, the present model has some 
similarities with the one-dimensional facilitated Ising model at zero temperature. Although the possible elementary 
dynamical events are different, both systems get eventually stuck in a metastable state characterized by all the holes 
being isolated, i. e., surrounded by two particles. 

The particle-hole description allows us to solve analytically the time evolution of the system. A closed hierarchy of 
equations can be written for the probability distribution functions D^. of finding r + 1 consecutive holes in the system. 
The general solution of this hierarchy has been derived. In the long time limit, the system reaches a state in which 
Dr = for all r > 1, i. e., all the holes are isolated. Moreover, the asymptotic value of the density of holes is not 
unique, but it depends on the initial configuration. 

Another similarity between the present model and the ISFM shows up in the context of granular media. When both 
models are applied to study tapping processes, the steady state reached by the system is consistent with Edward's 
thermodynamic theory [^ , p3| . This has been shown numerically in Lefevre and Dean's model ||2^ , and analytically 
for the ISFM In the latter, an effective dynamics for the tapping process has been derived. This allows to find 
explicitly the steady distribution function of the system, as well as the relationship between the vibration intensity 
and Edwards' compactivity. Therefore, it seems worth trying to derive the corresponding effective dynamics for the 
tapping process in the model considered in this paper 134]. 
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